Charge redistribution at YBa 2 Cu307-metal interfaces 
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Charge redistribution at interfaces is crucial for electronic applications of high-T c superconduc- 
tors, since the band structure is modified on a local scale. We address the normal-state electronic 
structure of YBa2Cu307 (YBCO) at an YBCO-metal contact by first principles calculations for 
prototypical interface configurations. We derive quantitative results for the intrinsic doping of the 
superconducting CuC>2 planes due to the metal contact. Our findings can be explained in terms of 
a band-bending mechanism, complemented by local screening effects. We determine a net charge 
transfer of 0.09 to 0.13 electrons in favour of the intraplane Cu sites, depending on the interface 
orientation. 
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Electronic transport in wires and tapes from high-T c 
materials is seriously affected by structural defects and 
interfaces. In contrast to most conventional superconduc- 
tors, bending of the band structure due to local variations 
of the charge distribution is strong enough to control the 
transport properties [l|, 0] • This results from large dielec- 
tric constants and small carrier densities, characteristic 
for high-T c materials The Thomas-Fermi screening 
length, over which band-bending is effective, therefore 
reaches the order of magnitude of the superconducting 
coherence length. 

The specific contact resistivity of YBCO-metal thin 
films is known to depend on the details of the contact ge- 
ometry In particular, the transport in micron-sized 
YBCO-metal hetero junctions is strongly affected by the 
orientation of the YBCO crystallographic axes with re- 
spect to the direction of the current flow [5|. From the 
theoretical point of view, effects of charge modulation at 
the surface of high-T c superconductors have been studied 
by Emig et al. [g]. It turns out, that surfaces are cov- 
ered by dipole layers, due to a local suppression of the 
gap function. Nikolic et al. j} study the charge imbal- 
ance at the boundary between a short coherence length 
superconductor and a normal metal by means of a self- 
consistent microscopic approach. However, first princi- 
ples electronic structure calculations taking into account 
the details of the crystal structure are missing so far, 
probably due to a high demand on CPU time. 

The technical optimization of interfaces calls for in- 
sight into the electronic structure close to the contact, 
which we address in the following for characteristical 
YBCO-metal interfaces. Since the electronic properties 
depend on the local atomic configuration, it is necessary 
to start from the details of the crystal structure in order 
to obtain reliable results. We present findings of band 
structure calculations for two prototypical YBCO-metal 
interface configurations fulfilling this requirement. Our 
calculations are based on density functional theory and 
the generalized gradient approximation, as implemented 
in the WIEN2fc program package [8|. This full-potential 
linearized augmented-plane-wave code is known to be 



particularly suitable for dealing with structural relax- 
ation and charge redistribution in complex geometrical 
arrangements [l(| HH ■ We obtain quantitative results 
for the band-bending magnitude and, therefore, for the 
intrinsic doping of the superconducting Cu02 planes due 
to the metal contact. 

Since band-bending is proposed to take place on the 
length scale of the YBCO lattice constant, the electronic 
structure of YBCO-metal interfaces becomes accessible 
to a supercell approach with periodic boundary condi- 
tions. In the following, we address supercells both par- 
allel and perpendicular to the crystallographical c-axis. 
In each case, we start from the experimental YBCO bulk 
lattice constants [12[ and optimize the atomic coordinates 
in order to minimize the atomic forces. In a second step, 
the structural relaxation of the supercells is carried out 
[HI, ■ Figure Q] shows the interface configurations un- 
der consideration, which we call the parallel (a) and the 
perpendicular (b) interface, referring to their orientation 
with respect to the Cu02 planes. 

It is convenient to choose fee Pd as the metallic sub- 
stituent, because of a minimal lattice mismatch of about 
0.7%. The supercell for the parallel interface consists of 2 
YBCO unit cells, whereas 3 YBCO unit cells are used for 
the perpendicular interface. The YBCO domain in each 
case terminates with a Cu-0 layer [l^, [1(1 [TtJ ■ Further- 
more, the metal domain comprises 3 and 4 Pd fee unit 
cells for the perpendicular and parallel configuration, re- 
spectively. Relative shifts between the YBCO and metal 
domain parallel to the interface have not be taken into 
consideration, since they do not affect our further conclu- 
sions [lbT | . Structure optimization results in a strong ten- 
dency towards Pd-0 bonding, whereas repulsion is found 
between Cu and Pd atoms. Importantly, the structural 
relaxation affects almost only the first atomic layer off 
the contact. Bond lengths calculated for these sites are 
summarized in Table [J for both supercells. 

In the following discussion, we compare the electronic 
structures obtained for our YBCO-metal interfaces with 
the results of a bulk YBCO calculation. We therefore 
mention that our bulk YBCO density of states (DOS) 
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FIG. 1: YBCO-metal interface configurations: (a) parallel 
and (b) perpendicular. The naming refers to the orientation 
of the interface with respect to the CuC>2 planes. 



TABLE I: Selected bond lengths at the parallel and perpen- 
dicular YBCO-metal interface, as resulting from the structure 
optimization. 
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FIG. 2: Partial Cu 3d densities of states for Cu02-plane sites 
in the vicinity of the parallel (a) and perpendicular (b) inter- 
face. The gross shape of the interface DOS resembles the bulk 
DOS when the latter is shifted by 0.20 eV (a) and 0.15 eV (b) 
to lower energies, respectively, which corresponds to charge 
carrier reductions of 0.13 and 0.09 holes per Cu site. 



agrees perfectly with previous theoretical and experimen- 
tal findings, see [3, 0, [2(| and the references given 
therein. Figure [5] shows partial Cu 3d DOS curves for 
Cu sites in the Cu02-planes of bulk YBCO and our in- 
terfaces (a) and (b). In the latter cases, the results refer 
to the Cu sites second closest to the contact. As com- 
pared to bulk YBCO, additional Cu states appear in the 
vicinity of the Fermi energy for the Cu atoms actually 
forming the interface, which trace back to a modified 
Cu-0 bonding. In contrast, the second closest Cu atoms 
resemble the bulk YBCO DOS very well and hence allow 
us to study the band-bending. For the parallel interface, 
we observe almost perfect agreement of the DOS curves, 
whereas some effects of the structural relaxation are left 
for the perpendicular interface. However, in both cases 
the bulk DOS has to be shifted to lower energies in or- 
der to reconcile the curves. The necessary shift amounts 



to 0.20 eV for the parallel interface, whereas 0.15 eV are 
sufficent in the perpendicular case. 

Since the gross structure of the Cu 3d DOS is not af- 
fected by the interface, we can interpret these energeti- 
cal shifts in terms of almost ideal down bending of the 
electronic bands due to a modified Fermi level. As a con- 
sequence, the hole count at the Cu sites is altered. To 
be specific, a shift of 0.20 eV corresponds to a reduction 
of 0.13 holes, and a shift of 0.15 eV comes along with a 
loss of 0.09 holes. Both these values and the calculated 
magnitudes of the band-bending are expected to be inde- 
pendent of the metallic substituent used for the interface, 
which we have confirmed for silver contacts. Importantly, 
the charge transfer likewise depends only little on the ori- 
entation of the YBCO-metal interface with respect to the 
unit cell of the high-T c compound. An intrinsic doping 
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close to 0.1 eV hence appears to be a general result for 
YBCO-metal contacts. Of course, oxygen defects close 
to the interface could modify the charge transfer. 

Core levels at atomic sites near the interface show en- 
ergetical shifts of about twice the magnitude reported for 
the Cu valence states. To understand this fact, electronic 
screening has to be taken into consideration. Screening 
is more efficient for the Cu 3d states than for any core 
states, since the former have finite weight at the Fermi 
energy. In general, the electrostatic screening length is 
only a few nanometers in high-T c cuprates and the crys- 
tal structure is very inhomogcncous. Conventional band- 
bending models based on a continuum description of the 
charge distribution thus cannot be applied. Nevertheless, 
electronic screening results in a significant reduction of 
the band-bending magnitude at YBCO-metal interfaces, 
therefore in a reduced charge transfer. 

Xu and Ekin report on specific resistivities for YBCO- 
Au interfaces of 10 -4 fi cm 2 to 10 -3 cm 2 at low temper- 
atures [2lJ. However, our calculations for normal-state 
YBCO-metal interfaces do not show a relevant reduction 
of the Cu 3d DOS at the Fermi energy. Even though no 
insulating layer is formed in the vicinity of the YBCO- 
metal contact, the observed interface resistivity can be 
explained in terms of screened band-bending. Since the 
charge carrier density in the Cu02-planes is significantly 
reduced, a local breakdown of the superconductivity is 
to be expected. 

We have presented electronic structure calculations for 
prototypical contact configurations between the short co- 



herence length superconductor YBa2Cu3 07 and a normal 
metal. In particular, we have discussed the charge redis- 
tribution in high-T c materials induced by normal metal 
interfaces. Calculations for well-relaxed supercells show 
that the charge redistribution can be interpreted in terms 
of an intrinsic doping of the superconductor on a nanome- 
ter length scale, in correspondence with the experimental 
observation of charge carrier depletion. The net charge 
transfer in favour of the copper sites in the Cu02-planes 
amounts to 0.13 electrons when the interface is oriented 
parallel to the superconducting planes, and to 0.09 elec- 
trons for the perpendicular orientation. 

This weak dependence of the charge transfer magni- 
tude on the orientation of the interface lets expect that 
the intrinsic doping at YBCO-metal interfaces is given by 
some 0.1 electrons for any contact geometry, as long as 
the YBCO domain terminates with a Cu-0 layer. More- 
over, the mechanism of screened band-bending and the 
net charge transfer are almost independent of the specific 
high-T c material and metal forming the contact. There- 
fore, the results are very general and can be applied to a 
large variety of interfaces. 
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